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Using an effective potential method, a replica formulism is set up for describing supercooled liquids 
near their glass transition. The resulting potential is equivalent to that for an Ising spin glass in a 
magnetic field. Results taken from the droplet picture of spin glasses are then used to provide an 
explanation of the main features of fragile glasses. 
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Structural glasses have been studied extensively at 
mean-field level by relating them to infinite-range p- 
spin spin glass models 0. The dynamical equations of 
these spin glass models are similar to those of the mode- 
coupling theory of liquids [2J • These equations predict a 
dynamical transition at a temperature Td below which an 
ergodicity breaking occurs. In the infinite-range p-spin 
glass models, a discontinuous thermodynamic transition 
also takes place at a temperature Tk and is usually iden- 
tified with the Kauzmann glass transition 0. T K < T d . 
It is expected that outside the mean-field limit, there 
will be no genuine dynamical transition at Td due to the 
fact that the metastable states which trap the system 
dynamically are unstable in finite dimensions because of 
activation processes over the finite free-energy barriers 
which separate the metastable states. However, it is com- 
monly assumed that the thermodynamic glass transition 
Tk could still exist in principle in finite dimensional sys- 
tems. This paper is essentially an investigation of this 
question. 

We start by constructing an effective Hamiltonian H[q] 
for supercooled liquids near the putative thermodynamic 
glass transition in finite dimensions. We find that H[q] is 
given by the same generic functional which describes the 
one-step replica symmetry breaking transition for the 
p-spin spin glass model. Our derivation is based on the 
effective potential method developed by Franz and Parisi 
for spin glass models. A similar effective Hamiltonian 
to ours was used in Refs. and for structural glasses, 
but not explicitly derived. 

The next step is the elucidation of the properties of 
H[q\. There are two coupling constants in the theory, wi 
and W2 and at mean-field level, the thermodynamic tran- 
sition is discontinuous (first order) if their ratio wijw\ 
is greater than one and continuous when less than one. 
There are arguments for a growing lengthscale in struc- 
tural glasses as the temperature is reduced (0,0, Ref. 
and references therein, and Bouchaud and Biroli 0) 
and experimental evidence for a lengthscale of several in- 
terparticle diameters already exists |9( . We feel that it is 
more natural to study parameters appropriate to the con- 
tinuous transition rather than the first-order transition in 



order to explain a growing lengthscale. Most studies of 
p-spin models have focussed on parameters which gener- 
ate the first-order transition to the one-step replica sym- 
metry breaking state, simply because in that situation 
one also has a mechanism for explaining the dynamical 
transition Td- However, it is clear from our derivation of 
the p-spin functional that the q a /3 in it is not simply re- 
lated to density fluctuations so our functional cannot be 
used to directly study dynamical effects in density corre- 
lations functions. In our view, they are best determined 
by mode-coupling theory 0. Our p-spin functional is 
appropriate only for the study of the putative glass tran- 
sition at Tk- 

Numerical studies have been performed on a p-spin 
model with p = 3 by Parisi et al. j_which at mean- 
field level has a continuous transition [ll|, (because for 
it the ratio w^/vji < I). They found a rapidly grow- 
ing lengthscale in this p-spin model as the temperature 
was reduced and an attempt was even made to obtain 
its critical exponents in four dimensions. Moore and 
Drossel 12j later showed that the continuous transition, 
if any, would be in the same universality class as the Ising 
spin glass in a field - the transition usually referred to 
as the de Almeida- Thouless (AT) transition [l^. How- 
ever, it has been argued that the AT transition does not 
actually exists for systems whose dimensionality is less 
than six [J^j, but instead there is a crossover line in the 
field-temperature phase diagram where the growing cor- 
relations between the spins saturate. The value of this 
lengthscale can be estimated from droplet scaling argu- 
ments 15]. In the numerical study of Parisi et al. |l(J 
the largest system size studied was of linear dimension 
six lattice spacings and if the correlation length is larger 
than the size of the system then the crossover would look 
like a phase transition. In the work of Moore and Drossel 
[lij ] the lengthscale in the (three-dimensional) Migdal- 
Kadanoff approximation saturated at around 28 lattice 
spacings at low temperatures. We shall also show using 
droplet ideas borrowed from spin glasses how this grow- 
ing lengthscale can produce dynamical behavior of the 
Vogel-Fulcher type [lf| . 

The calculations start with defining an overlap p c (r) = 



2 



5pi(r)Sp2(r) between two configurations of density fluc- 
tuations Sp = p — (p) in two copies of the liquid. By av- 
eraging over the density configurations of the first copy 
while holding the density configurations in the second 
copy fixed at fixed values of the overlap p c (v) we obtain 
a constrained partition function 

Z[ Pc (r),Sp 2 (r)} = (S(p c (r)-Sp 1 (r)Sp 2 (r))) Pl . (1) 

The effective potential is given by the average of the free 
energy with respect to the density configurations in the 
second copy: S![p c (r)] = — T(ln Z[p Cl Sp 2 ]) P2 . In order to 
perform the average of the logarithm, we use the replica 
method and write \nZ = lim n ^o(Z n — l)/n. Using the 
integral representation of the delta function in l|TJl. we 
can write the effective potential as the following averages 
over p 2 and pf (a = 1, • • • , n): 

il\p c (T)]=-T /il^ ex P *E / drA -(r)Pc(r 

x(/exp -i^T J drSp"(r)Sp 2 (r)X a (r 

We evaluate the above averages over pf and p 2 by 
making successive cumulant expansions. The result- 
ing integral can be expressed in terms of the liquid- 
state correlation functions G(ri,r2) = (6p(rx)6p(r2)), 
G(ri,r 2 ,r 3 ) = {5p{Y 1 )5p{Y 2 )5p{rz)) 1 etc. Keeping terms 
up to 0(A 3 ), we have 



Q[p c (r)] ~ -T /ll^ ex P *'E / drX *(r) Pc (r 
~ f dld2G 2 (l,2)^A Q (l)A Q (2) 
~ J dld2d3 G 2 (1,2,3)^A Q (1)A Q (2)A Q (3) 




x exp 



x exp 
xexp [0(A 4 Q ,A 2 A 2 )], 



(3) 



where 1,2, ■ • • denote the position vectors ri, r 2 , ■ ■ ■ . 

Now we define q a /3( r ) = A a (r)A/3(r) for a ^ [3 and 
write the effective potential as a functional integral of 
the fields q a p. We insert the following identity inside the 
above expression: 



-/iWn^r 

J a<(3 J a<f3 



(4) 



x exp 



dr u af3 {v) (q a p(r) - A Q (r)A /3 (r)) 



Due to the coupling between A Q and u a 0, the integration 
over X a generates various terms in u a p- For example, to 
the lowest order, the quadratic term in u a p comes from 



joining two pairs of A's in the two wAA vertices. Therefore 
the coefficient of the quadratic term in u a p involves, to 
the lowest order, the square of the A-propagator K(l, 2) 
satisfying 



j d3K(l,3) [G 2 (3,2)] =5(12). 



(5) 



There are, of course, other contributions to the 0(u 2 ) 
term, which come from the higher order vertices in A. 
These just renormalize its coefficient. Two kinds of cu- 
bic terms in u a p are generated, namely, Tr u 3 and u c\n- 
The coefficients of these terms will contain in general the 
contributions from the higher order vertices in A. Keep- 
ing up to the cubic order in u a p, we have 



nMrw ~ t [ n v qa , / n 

xexp[i^ / dr u a p(r)q a p(r)] 

a<f3 J 

x ex P [\ Yl I dld2d3 2 ' 3)w Q/3 (l)Pc(2)p c (3)] 
x exp I — - 

/ dld2 B(l,2)u a0 (l)u af3 (2)] 

a<f3 J 

xexp[^ dld2d3 G(l, 2, 3)u Q/3 (l)u^ 7 (2)u 7Q (3)] 

(«,/3,7) 

x exp 

[\zZj dld2«H-D(l,2,3)« o/3 (l)u a/ j(2)u a/ ,(3)], 



(G) 



a</3- 



where we dropped the terms in p c which do not couple 
to u a p. The coefficients A, B, C and D are given, to the 
lowest order, by 



4(1,2,3) ~A-(1,3)JC(2,3), 

S(1,2)~A- 2 (1,2), 

G(l,2,3)~X(l,2)tf(2,3)X(3,l), 



(7) 
(8) 
(9) 
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L>(1,2,3)~- / ndr 2 G 2 (4,5,6)G 2 (7,8,9) (10) 

^ i=4 

x K(l, 4)K(1, 7)K(2, 5)K(2, 8)K(3, 6)K(3, 9). 

Finally we integrate over u a p to obtain the effective 
Hamiltonian in terms of q a p as follows. First note that 
in addition to the usual quadratic term in q a p, there 
are cubic terms in q in the form Tr q 3 and g 3 ^ with 
real coefficients. We also note that the terms involv- 
ing p c couple to q a p in the form p^q^a with integers 
k and k' . The resulting expression gives the effective 
potential f2[p c ] in terms of a functional integral over 
the fields q a p. Since the value of the overlap is deter- 
mined from the condition 5£l/5p c = 0, p c (r) = is al- 
ways a solution and describes the liquid phase. We take 
the effective Hamiltonian H[q] for supercooled liquids as 



3 



il[p c — 0] ~ fY\ a< pT>q al sexp[—H[q]]. The coefficients 
of various terms in q a p in H[q] are in general expressed 
in a nonlocal way much like Eqs. 0-JTOJ. However, 
the most relevant terms will be the local ones where all 
the points in space are the same. Keeping as usual the 
derivative term only in the quadratic order in q a p, we 
obtain the effective Hamiltonian as 

a</3 a<f3 

-!^(,)-!££^(r)}. (n) 

a</3 

The coefficients c, t, w\ and W2 will be functions of the 
temperature and density of the liquid, with smooth de- 
pendence on them. The study of the first-order transition 
which occurs at mean-field level requires the inclusion of 
terms quartic in the q a p fields, but these are 'irrelevant' 
variables at the presumed continuous transition in finite 
dimensions and they will be dropped. At mean-field level, 
the transition is continuous if 1112/ wi < 1 [4|. In the p- 
spin replicated functional of Eq. l|llfl the coefficients in 
it have contributions from all the higher order cumulants 
of the density etc.. Fortunately as we are studying a con- 
tinuous transition, their numerical values will only effect 
non-universal features of the putative transition. 

The thermal average {q a j3) is non-zero even above any 
supposed transition (when one goes beyond the mean- 
field approximation) because of the term in W2 in the 
Hamiltonian of Eq. As a consequence the only pos- 

sible continuous phase transition is one which breaks the 
replica symmetry of the high-temperature phase. This 
transition was shown in Ref. [12j to be in the same uni- 
versality class as that for the transition from the param- 
agnet to the spin glass phase in a field, first discussed by 
de Almeida and Thouless [l3j. The Hamiltonian of an 
Ising spin glass in a field is 

H = - Jij SiSj ~hJ2s l . (12) 

<ij> i 

The couplings Jy are quenched random variables con- 
necting (usually) nearest-neighbor sites. The Edwards- 
Anderson order parameter is q — ^2i(Si) 2 /N and in zero 
field becomes non-zero below T s , the critical tempera- 
ture of the zero-field spin glass transition. Using the 
replica method one can set up a functional analogous to 
that for structural glasses in Eq. i|ll|) for the spin glass 
Hamiltonian [T^J. It differs only in that the cubic term 
(W2/3) ^2 c?„0 is absent but there is instead a linear term 
h? ^2 Qa/3- Now when the temperature T is in the vicin- 
ity of the thermodynamic structural 'glass transition' Tr-, 
we expect it to be well below T s . At such temperatures 
W2 and h 2 can be related by noting that the cubic term 
could be approximated by w 2 {^ptlap- (lap) wul be 
independent of a and j3 as we are working in the replica 



symmetric phase, and only slowly varying as a function 
of temperature so then h 2 cx W2 ■ Then because W2 is a 
function of the temperature and density of the liquid, h 2 
in the Ising spin glass analogue will also be a temperature 
and density-dependent quantity. 

The question then of whether there is or is not a struc- 
tural glass transition then turns on whether there is an 
AT line in the spin glass analogue. This has been con- 
troversial in the past, but numerical [isj | and experimen- 
tal data 0| favors its absence in three dimensional spin 
glasses and we further argued in Ref. 0] that for all di- 
mensions d < 6 there would be no AT line. We shall 
henceforth take it that there is no AT line and hence 
that there is no thermodynamic glass transition for three 
dimensional systems. How then can we explain the phe- 
nomena normally associated with the glass transition? 

The key to understanding this point is to realize that 
even if there is no AT line there can be a growing length- 
scale as the temperature is reduced. The lengthscale is 
analogous to the 'mosaic' lengthscale of other authors 
Q El lL S> but its physical origin is different. In our 
case, the lengthscale is equal to the size of the droplet 
in the spin glass analogue which is typically flipped by 
the application of the field h. To flip a region of size 
L in the low-temperature phase of a spin glass typically 
costs an 'interface' free energy of order (L / £,(T)) e , where 
£(T) is the correlation length associated with the spin 
glass transition in zero-field and 9 is an exponent typi- 
cally ~ 0.2 15). Each droplet has a magnetic moment of 
order \J~qL d l 2 so in a field of magnitude h the domains 
are of size L*, where 

(L*/aT)f ~ V^(L*) d/2 h, (13) 

in order that the energy gained by flipping a droplet in 
the field h just balances the interface free energy cost. 
This result for L* should be valid provided L* > £(T). 

In the droplet picture the free-energy barrier which 
has to be surmounted in order to flip a region of size 
L* is A(L*/£(T))^, where A is of order of the critical 
temperature T s |l5|. The basic relaxation time r in the 
system is of Arrhenius form: 

T~T exp[A(L*/Z(T)^/T}. (14) 

Alas, the value of ip in three dimensional spin glasses 
is very poorly determined, with quoted values varying 
widely. 

When T <C T s , £(T) will be a slowly varying function 
of temperature and of order of the interparticle spacing. 
q will also then be nearly temperature- independent. If 
W2 is only weakly temperature dependent, L* will be ap- 
proximately temperature independent and r will then be 
of the simple Arrhenius form. This will be the case then 
of strong glasses. For fragile glasses we expect that 

w 2 ~ h 2 ~ (T - T ), (15) 



4 



so that ui2 appears to be going to zero at T . We stress 
that the effective vaiue of wi probably will never van- 
ish, since even if the bare coefficient did go to zero the 
higher terms in the functional H[q] which break the time- 
reversal invariance will make its effective value non-zero 
so Eq. H15fl i s j us t a simple phenomenological way to 
obtain a growing lengthscale. Then 



tion of this paper. We are indebted to Drs G. Biroli and 
J-P Bouchaud for their comments and insights. 



T exp 



C 



(16) 



where C denotes a constant, ft is clear that this has 
similarities to the Vogel-Fulcher 0] form for structural 
relaxation times in the supercooled liquid: 



r exp 



C 



T —Tn 



(17) 



Note that we have proceeded on the assumption that the 
exponent Tp for the supercooled liquid and the spin glass 
are the same. In principle our mapping between liquids 
and spin glasses is explicit only for 'equilibrium' quanti- 
ties and might not extend to dynamical features. Note 
that the droplet formula for r of Eq. (|H|) is only ex- 
pected to hold when L* is much larger than the typical 
distance between the molecules of the glass, and exper- 
iments might not have reached the region where it has 
become valid; in other words the simple exponent of the 
Vogel-Fulcher fit might reflect pre- as ymp totic behavior. 

The simulations of Parisi et al. [l0( show that the 
entropy per spin went to zero apparently as s c (T) « 
fcs(T — Tk), although no data was obtained close to 
Tk and without a genuine transition the entropy is un- 
likely to ever vanish completely. A similar behavior is 
found for the configurational entropy per particle in frag- 
ile glasses. What is striking about the simulations is that 
the temperature Tk was very close to the temperature To 
at which the correlation length was apparently diverging 
(the difference was less than 1%). One feature common 
to fragile glasses is that for most of them To extracted 
from the Vogel-Fulcher fit is close to the Kauzmann tem- 
perature Tk, the temperature at which the extrapolated 
entropy of the supercooled liquid apparently falls to that 
of the crystalline phase. The ratio Tk/Tq lies between 
0.9-1.1 for many glass formers amongst which Tk ranges 
from 50 K to 1000 K j2jj. Since we are arguing that 
there is no genuine transition but only a crossover, the 
'T ' obtained from the timescale t, the T of Eq. lfl"5|l 
and Tk need not be precisely the same. 
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